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Laboratory for Planetary Atmospheres 
ABSTRACT 

The characteristic enhancements of ring current particles 
with energies of about 1 to 100 keV, associated v/ith mag- 
netospheric substorms, were observed by Explorer 45 
(S^-A) around the plasmapause in the afternoon to midnight 
region, and showed the v^haracteristic structure called a 
'nose' in the proton spectrograms. This paper describes the 
motion of these particles in the equatorial magnetosphere, 
under a recently proposed convection electric field and a 
dipole magnetic field. Approximate equations of a bounce 
period, a second adiabatic invariant and a bounce-average 
azimuthal velocity are given with inaccuracies less than about 
10“^ for all pitch angles. The complete set of flow patterns 
of 90° pitch angle particles is also presented by means of 
stagnation lines through which radial drifts and/or azimu- 
thal drifts change their directions . The particle tracings in 
the magnetosphere give a basic concept to explain the ob- 
served nose characteristics. 


i i i 


CONTENTS 


U 

\\ 


Page 

ABSTRACT iii 

I. INTRODUCTION ^ 1 

II. A PREMISE 2 

III. MOTION OF PARTICLES WITH ARBITRARY PITCH ANGLE . 4 

Pitch Angle Change . . . 6 

Particle Energy Change 7 

Approximations for f(y), I(y) and G(y) 7 

IV. ELECTRIC FIELD MODEL AND PARTICLE TRAJECTORIES . 9 

V. MOTION OF 90° PITCH ANGLE PARTICLES 12 

Stagnation Lines and Last Closed Equipotential Line 13 

Particle Flow Patterns With 7 -- 2 . 16 

VI. TIME-DEPENDENT PARTICLE FLOWS 20 

VII. ACKNOWLEDGMENTS 22 

VIII. REFERENCES 23 

APPENDIX A . . 29 

APPENDIX B 31 

APPENDIX C 33 


PAGE Fima, 


V 



TRAJECTORY TRACES OF CHARGED PARTICLES 
IN THE MAGNETOSPHERE 


I. INTRODUCTION 

There have been a number of observations of charged particles in the magneto- 
sphere related to magnetospheric substorms [e. g. , Vasyliunas, 1968; Shield and 
Frank, 1970; Lezniak and Winckler, 1970; DeForest and McHwain, 1971; Frank, 
1971], It has been recognized that the convection electric field plays an impor- 
tant role in the motion of those particles [e.g. , Taylor and Hones, 1965; Nishida, 
1966; Kavanagh et al. , 1968; Mcllwain, 1972; Gumett and Frank, 1973; Heikkila, 
1974; Smith and Hoffman, 1974] . Mcllwain (1973) deduced his electric field 
model with a specified injection boundary [Mauk and McHwain, 1974] to explain 
ATS-5 particle data associated with substorms. Konradi et al (1975) followed 
McHwain’s model and showed thgt the substorm associated particles observed by 
Explorer 45 came from along or beyond the Injection boundary. He also calcu- 
lated the 27° pitch angle proton dispersion pattern and compared it with the 81° 
pitch angle protons. 

Many theories have been developed about the charged particle motions in the mag- 
netosphere under the convection electric field. Chen (1970) gave a morphology 
of 90° particle flow patterns assuming a uniform dawn-dusk electric field (see 
also references in his paper for other works). Stem (1975) calculated the par- 
ticle energy change and showed some particle trajectories under the recently 
proposed convection electric field [ Volland, 1973; Stem, 1974], As for the off- 
90° pitch angle particle motions at the equator, several approximation methods 
to compute particle bounce period, second invariant and bounce average drift vel- 
ocity have been investigated [Lew, 1961, Hamlin et al. , 1961; Lencheck et al. , 
1961; Schultz, 1971; Chen and Stem, 1975], But the off-90° particle trajectories 
in the magnetospheric equatorial plane have not been intensively studied yet. 

Smith and Hoffman (1974) observed with Explorer 45 (S^~A) the storm-time ring 
current proton enhancements and found in the 90° pitch angle proton spectrograms 
the characteristic ’nose’ stmctures inside the plasmapause. They concluded that 
those particles came to the lower L-value region from the plasma sheet or tail 
regions under convection, corotation, and magnetic gradient and curvature drifts. 
An example of such data appears in Figure 1 in the form of spectrograms at eight 
different pitch angles for both outbound and Inbound passes during the magnetic 
storm of February 13, 1972. (For a description of Explorer 45 and the particle 
detectors, see Williams et al. , 1968, and Longanecker and Hoffman, 1973). 

Since the orbital plane was almost equatorial, pitch angles averaged over about 10 
degrees. From the DC electric field experiment [Maynard and Cauffman, 1973; 
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Cauffman and Maynard, 1974], the plasmapause positions during the outbound 
and the inbomid passes were identified as L = 5. 1 and L = 4. 8, respectively. It 
is obvious from these data that the typical ’nose' structure in each pitch angle 
spectrogram was observed inside the plasmapause. And, these structures def- 
initely depend on particle energy and their pitch angles and the time-phase in 
the substorm. The higher pitch angle particles were detected at the lower L- 
value positions along both outbound and inbound orbits. The energy of the 'nose' 
edge was about 25keV in the outbound and about 13 keV in the inbound for all pitch 
angles. Before the enhancements of energetic particles were observed, there 
were several sub storm onsets identified from both high-latitude and low-latitude 
ground magnetograms. 

To interprets these particle enhancements, especially the lkeV-50keV proton 
populations associated with substorms, we speculate the particle sources in the 
tail region which start to move towards the earth due to an enhancement or mod- 
ification in the convection electric field. The simultaneous differential equations 
are deduced to present the motion of these particles in the equatorial plane, and 
approximate formulas for bounce period, second adiabatic invariant and bounce- 
average azimuth drift velocity of arbitrary pitch angle particle are given with 
sufficient accuracies of about 10~^. Also, the 90° pitch angle flow patterns are 
delineated by means of stagnation lines (the definition of a stagnation line is given 
in V) through which particle flow directions are changed, giving complete mor- 
phology of particle motions in the equatorial plane by using the recently proposed 
convection electric field model of Volland and Stern [Volland, 1973; Stem, 1914, 
and 1975; Maynard and Chen, 1975] . The examples of particle tracings including 
off-90° pitch angle particles in the magnetosphere give a basic concept to explain 
the newly injected particle populations in the magnetosphere. 

MKS units are adopted through this paper unless otherwise mentioned, and the 
symbols for quantities relevant to charged particle motion are listed in Appendix 

A. 


II. A PREMISE 


The equation of charged particle motion in a magnetic field B and an electric 
field E is given by 



q 
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where r is the position of the particle, q the^ particle charge, m the relativistic 
mass, and F' non-electromagnetic forces such as gravitational force, inertial 
force, gradient of pressure etc. Here we develop the above equation in a region 
where external forces F’ may be neglected and a particle motion associated with 
a substorm injection may be well described by the non-relativistic equation (m = 
my, tlie rest mass of the particle). Furthermore, we assume the following 
conditions; 

1. the geomagnetic field is a dipole field, its axis being perpendicular to 
the equatorial plane, and a field line is an equipotential, 

2. changes in E and B are very slow, and the particle motion can be rep- 
resented as a motion of its guiding center (cyclotron average drift 
motion, e.g. , Northrop, 1963), 

3. the particle motion is adiabatic, i. e. , first and second adiabatic invari- 
ants iiJi and J) are conserved, 

4. there are no local energization and loss processes. 

From the above basic equ'atlon, a guiding center drift velocitj^ iTj, in the equatorial 
plane is obtained, by ave 7 jaging over a cyclotron motion and a bouncing motion 
between particle mim rcx' points. 


iTp = Fj^xB/qB^ (1) 

and 

Fj^ = q • E - q(co x R) x B - W • GCUq ) • Vj^B/B, (2) 

where R is the position vector from the center of the earth in the equatorial 
plane, W the particle kinetic energy, cj the angular velocity of the earth, the 
particle pitch angle at the equator, and the function G(ao) is given later. 

The first term is a force due to an external electric field, the second term rep- 
resents the corotation of the particle with a magnetic field line, and the tliird 
term shows the combined magnetic gradient and curvature drift. In the case 
a ^ 0°, the particle kinetic energy W is written as 

W = mB/y“ (3) 

where y = sinup, and \x = Wj^/B = l/2mv^/B, the first adiabatic invariant. If 
(Xg = 90°, the third term becomes a simple form. 



Note that the conservation of the first invariant is of no avail if = 0“ (y = 0). 

It is clear from equation (1) that the mndamental particle velocity vector direc- 
tion in the midnight region is (i) earth-ward under a dawn-dusk convection elec- 
tric field, (ii) eastward under the corotation field, and (iii) westward for protons 
and eastward for electrons under gradient- curvature drift. Combinations of 
these three kinds of drift motions make complicated particle flow lines which def- 
initely depend on particle energy and pitch angle, 

m. MOTION OF PARTICLES WITH ARBITRARY PITCH ANGLE 

To solve equation (1) with initial conditions (R-, y^ Wj; t = 0) in given electric and 
magnetic fields, E = E(R, P) and B = B(R, P), where P is magnetic local time, 
we have to obtain the kinetic energy W or pitch angle y (= sin ) as a function of 
a particle position R(R, P), and a functional form of G(y) which is closely re- 
lated to a particle bouncing motion in a mirror field. 

There are three fundamental equations in connection with particle motion in a 
dipole field [e. g, , Hamlin et al. , 1961; Roederer, 1970] , as follows: 

1, Bounce period ; 




2 ds 

Os' 

m 


4R 

V 




(5) 


where 


and 


f(«o) 



cos X(1 + 

(l -siii^aQ • sec^X(l + Ssin^X)'"^^} 


dX, 


( 6 ) 


cos'^X^ 


2 — m 


(I + SsiiOX^)’/^ 


(7) 


2. Second adiabatic invariant in the case of no parallel external force; 


J = <^p„ ds = 2pRl(o:„) 


(8) 



where 



sin-ttg • sec^Xd + (1 +3sin^xy^^ cosXdX. 


(9) 


3. Bounce-average azimuthal drift velocity; 


d(j) _ Sinv^ ^ g(a„) 

dt 2qk'o f(«o) ’ 


(10) 


where 


g(0!o) 


r • sec^XCl + 3sin-X)^^0 (1 +simX) cos^X 

Jo n + 3sin^X)^/2 _ s{n2Q,^ . sec6x(] + Ssin^X)'/^} 


(11) 


and 

kp = a^kg = B/R3. 


After manipulation of equation (10), the third term of equation (2) can be obtained 
andG(ao) = g(ao)/f(o;o). 


To understand the particle motion we must first compute an equatorial pitch angle 
change and a particle energy change due to radial drift motion. From equations 
(6) and (9), there is a relation of f(ao) I(c>;o) ^ differential form as 


dy 




(12) 


where 


y = sinay. 

On the other hand, the constant of ju/j^ gives 

!(:. ) L.^ 

^ ' l(yj) L’/^ 

Suffix (i) denotes initial values. 


( 13 ) 
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And, the conservation of ju gives (when M 0) 


> 


fj 


” (Wy^L^) = 0. (14) 

dL 

The other form of energy change is given by Northrop and Teller (1960), and 
Kivelson and Southwood (1975), 

dW/dR = qRB(d0/dt). 


This equation can be transformed to a simple form as 


dL 


W 

-3G(o;„) - ■ 


(15) 


Equations (14) or (15) give the energization rate of particle due to radial drift, 
which increases with a pitch angle increase. 


Pitch Angle Change 


For the equatorial pitch angle, Southwood and Kivelson (1975) suggested the 
approximate expression of 

y ^ yi 

This is used as an initial condition for solving equation (13) by the Wegstein iter- 
ation method, alter slight modification so that the right hand side of equation (13) 
never exceeds 1, i.e. y = sina^ < 1. This method has its advantage that the 
desired accuracy can be easily specified, if the value of I(y) is obtained with the 
same accuracy as discussed later. Figures 2a and 2b show changes in equatorial 
pitch angles, initial pitch angle being set at L = 10. The particles change their 
pitch angles less than 4 degrees in the region of L = 5 to 3. The calculation of 
particle trajectory, therefore, is done by tracing the same particle in spite of 
its pitch angle change as is shown in the figrires, since the experimental data in 
the case of Explorer 45, for example, are averaged over the pitch angle interval 
of about 10 degrees (see Fig. 1). To compute a particle velocity vector, however, 
the change in pitch angle must be taken into account. Note that an equatorial 
pitch angle o!q must satisfy the condition sinaQ > sinag = 1/ {L^(4 - S/D^^} 
where «g is an equatorial loss cone if we assume the mirror points above thq 
earth’s surface. This loss cone limit is also shown in the figures. 
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Particle Energy Change 


The particle energization due to radial drift motion is calculated using equa- 
tion (14) combined with equation (13) and the result is illustrated in Figrire 3. 
The limiting cases for = 90° and = 0° are given simply as energy changes 
of Wj^ and W„ from first and second adiabatic invariants, respectively: 


r 


Wf/Wii = (Li/L)3 for = 90°, 
W||/W| 1 j = (Lj/L)^ for = 0°. 


(16) 


The energy gain due to radial drift of the particle started from L = 10 is about 
10 times in the region around L = 4, and the difference in energy between 0° and 
90° pitch angle particles is 2. 5 times. As is obvious from the form of equa- 
tion (15), the energization rate is independent of initial energy. It is worthwhile 
to notice that, as for the gradient-curvature drift velocities of same energy 
particles with different pitch angles, the velocity of lower pitch angle particles 
is higher than the velocity of higher pitch angle particles in the region of high 
L- value. For example, in the case W(90°) = W(0°) at L = L’, W(0°)/W(90°) = 
(Lj/L') at L = Lj and then the third term of equation (1) or equation (10) gives 
#/dtloo/d<^/dt| 5 QO = (2/3) (Lj/L), since G(0°) = 2/3 and G(90°) = 1 (see Fig. 4). 
Therefore, \d<p(aQ = Q°)/dt | > | d0(o:Q = 90°)/dt| in the region where Lj>(3/2)L', 
in spite of ju(<^o ^ ^ " 90°). 


Chen and Stem (1975) gave an approximate formula of energy as a function of 
I(y)/y* Cowley and Ashour-Abdalla (1975), Nakada et al. , (1965), and others 
also show the changes in equatorial pitch angle and particle energy. To inter- 
prete the substorm associated particle enhancements in the ring current region, 
more important than the changes in pitch angle are the particle trajectories with 
time-development coupled with tlie changes in particle energy. 

Approximations for f(y), I(y), and G(y) 

In the calculations of particle trajectories by solving simultaneously the three 
differential equations as shown later, it is not practical to compute by integra- 
tion the function values of f(y), I(y), and g^) given by equation (6), (9), and 
(11), at every time step. 

There have been many approximation formulas obtained by Lew (1961), Hamlin 
et al. , (1961), Lencheck et al. , (1961), Schiilz (1971), and Chen and Stern (1975). 
Here a more accurate analytical approximation formula for f(y) is given, and 
I(y) and g(y) are deduced selfconsistently from f(y). Adopting Lencheck' s 
approximation formula of f(y) as a zero order equation and adding higher order 
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tex-ms of the form (n and m are integers), the best fit analytical expression 

is obtained in the case m = 3 and n = 1, 2, 3, and 4 in order to get the inac- 
curacy less than about 0. 1% for all pitch angles, i. e. , 

f(y) = a-jS(y + +a 3 Y + (17) 


where 


aj + a2 + as + a4 = 0 , 

a = 1 +2n(2+V3)/2V3 -1.38, 


(18) 


and 


|3 = (a/2 ) -(tt/1 2)7? ^0.32. 

This equation gives exact values at y = 0 (0° pitch angle) and y = 1 (90° pitch 
angle). To determine the values of the coefficients a, , . . , and a 4 by the best 
fit of the approximation function given by equation (17) to the function f(y) of 
equation (16), we must compute an accurate functional value of f(y), the inte- 
grand of which has a singularity, by developing a special numerical integration 
method (see Appendix B). The resxilts are; 

a, = 0.055,02 = - 0 . 037,03 - -0.074 and 04 = 0.056. (19) 

A similar approach has been used by Hilton (1971) to obtain Mcllwain''s function 
of L(I, Bjn) and by Chen and Stern (1975) for W, f, and g as a function of I(y)/ 5 ^ 

Using analytical expression of f(y) in equation (17), I(y) is easily calculated from 
equation (12) (Schulz, 1971): 

Ify) = 2 y r Lr 2 f(u)du 

Jy 

= 2a(l - y) + 2(3y2ny + 4|3(y - y'"^^) ( 20 ) ^ 

+ 3aj (y’ - y) + 6 a 2 (y^''^ -y) + 604 (y - 2 a 3 y 2 ny. 

From equations (12) and (13), we have the differential form of the pitch angle 
change with respect to time t, 

^ _ Z . 

dt ~ L * 4f(y) ’ 


c\L 

dt 


( 21 ) 
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Putting equations (21) and (15) into (14), the simple relation G(y), f(y), and I(y) 
is deduced as 


G(y) 



I(y) 
6f(y) ■ 


Using equations (17) and (20), g(y) can be written as 

3g(y) = 2a-(5(3-o:)y-(j3-a3)yCny-|3y’'^ 

+ I ajyi/3(i -yl/3) + 6a4y(yi/3 _i). 


( 22 ) 


(23) 


The computed functional forms of f(y), I(y), and G(y) are illustrated in Figures 
3(a-c) as functions of equatorial pitch angle and colatitude of mirror point. To 
discuss qualitatively die trajectory of particles with arbitrary pitch angle in the 
equatorial plane, these figures together with Figures 2 end 3 are very useful 
(see equation (1)), though an actual particle trajectory can be obtained with a 
sufficient accuracy using the approximate expressions of f(y), I(y), and G(y) in 
equations (17), (20), and (22). The accuracies of approximations of these func- 
tions depend on the equatorial pitch angle, and are shown in Figure 5; the rela- 
tive accuracy is calculated as 1 - (approximate value)/ (true value). All these 
approximate functions have inaccuracies of about 10"^, In the same figures are 
shown the results of Schulz’s approximation (1971) evaluated again using accurate 
values of f(y) and I(y). The great improvements of the new approximations 
are evident from the figures. The expressions derived by Chen and Stern 
(1975) have accuracies of about 1.7%, but their variable X = I(y)/y has a good 
advantage thatX* is constant, whereas our variable y should be calculated 
using equation (13) or equation (21). 


IV. ELECTRIC FIELD MODEL AND PARTICLE TRAJECTORIES 

Charged particle flow patterns in the magnetosphere were first analyzed inten- 
sively by Chen (1970), assuming a constant dawn-dusk electric field and a dipole 
magnetic field. Subsequently, the observed plasmapa.use and plasma trough 
locations by Ogo 4, Ogo 5, and Ariel 3 have been interpreted, introducing a 
time dependency of the electric field strength as a fimction of Kp (interplanetary 
magnetic activity index) [Chen and Greboy/sky, 1974; Grebowsky et al. , 1974; 
Chen et al. , 1975]. Mcllwain, on the other hand, proposed his own electric 
field model E3 (1972) and E3H (1973) including a corotation field due to a rota- 
tion of a dipole magnetic field. The observed particle data of ATS-5 associated 
with substorms can be explained in his electric field model and the distorted 
magnetic field model M2, 
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Recently, based on the observed electric field data (Heppner, 1972) and the shape 
of the plasmapause, Volland (1973) and Stern (1974, 1975) proposed the following 
large scale potential fields as 

= ART' sin 0. (24) 

Then, 

"E = -v4>£ (25) 

where A is a coefficient which determines the electric field intensity. Both of 
them reached the same exponent of 7 = 2. Maynard and Chen (1975) succeeded 
to interpret the locations of the plasmapause and isolated cold plasma regions 
observed by Explorer 45, by using this model with 7 = 2 , A being a function of 
Kp as 


A = 


0.045 

(1 - 0.159kp + 0.0093 kp2)3 


(kV/R2) 


(26) 


where Rj, is an earth's radius. They deduced this Kp dependence from the 
plasmapause observations by Ogo 3 and 5, and this produces cross polar cap 
potentials similar to those found by Heppner (1973). 

Using the corotation field expressed in equation (1), it can be readily obtained 
that the relation of a coefficient A in equation (24) with a stagnation point R = Rq 
( or L = Lq) at dusk for zero energy particles is 


and 


K 



K = -aTkQOj/Rg. 


(28) 


Hereafter, this electric field model represented in equations (25) to (28) is 
adopted. 


Introducing a dimensionless distance x (distance normalizedby Rq, i. e. x = R/Rq = 
L/Lg and Uj^ = Rq x), 'equation (1) can be written in polar coordinates (R, 0), 


d 

dt 


X 


dt 




+ x 


dt 




(29) 




f 
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where 


dx ^ 

“ = -CO -coscl), 

dt 7 


(30) 


and 


- — = x'>' + isin0 + 1 - 3W • G(y) 
CO dt ^ 




(31) 


The origin is the earth's center, and'^j^ and'^^ are unit vectors. Equations (30) 
and (31) together with equation (21) are fimdamental simultaneous differential 
equations of charged particle motion (the other forms of these equations are 
presented in Appendix C with conventional units). These equations do not con- 
tain explicitly the coefficient A of electric field strength which is only introduced 
in terms of scale factor Rq. From equations (27) and (28), 


Ro = 


'a^k^coy + 1 

. 7A / 


(32) 


In the case 7 = 2, 


Rn 

Lq = — ^ 3.6/(A[kV/R2])i/3_ (33) 

a 

The Kp dependencies of the stagnation point at the dusk and midnight crossing 
points of the last closed equipotential line of zero energy particles are shown in 
Figures, using equations (26), (32), and(54)in the case of7=2. In the steady elec- 
tric field model these are identical to the plasmapause locations. Mcllwain pro- 
posed a particle injection boundary from which energetic charged particles are in- 
jected into the magnetosphere at substorm onset time [Mcllwain, 1973; Mauk and 
Mcllwain, 1974; Konradi et al. , 1975]. The radial distance of this injection 
boundary Rj^ is given as R|j = (122-10 Kp)/ {(j) - 7.3) where Rj, is in L-value and 
<t> the local time in hours, and is illustrated also in Figure 6. This figure implies 
that Mcllwain' s injection boundary is far beyond the last closed equipotential line 
with an electric field strength used to explain the Explorer 45 plasmapause obser- 
vations by Maynard and Chen (1975). The shapes of the injection boundary and 
the last closed equipotential line in the dusk to midnight region is discussed later 
(see also Fig. 8). 


11 


V. MOTION OF 90° PITCH ANGLE PAETICLES 


From equations (30) and (31), fundamental equations of motion of 90° pitch angle 
particle (y = 1) are 


dx 

Tt 

1 ^ 

CO clt 


+ 2 

= -CO CQScj), 

1 


= ^ sin 0 + ] 


3ju' 


(34) 

(35) 


where 

X = R/Rg = L/Lg and - /u/qcoRg. (36) 

To find the particle flow lines tlae time is eliminated from equations (34) and (35), 
and by integrating with respect to x and 0, we have 


2 ^ 

7 


sin 0 + 




= c 


(37) 


where c is an integral constant. 

This equation, as a matter of course, is identical to the conservation law of total 
energy of particle; qd’ + [JB - C and C(=/cqc) is an initial total energy conserved 
as a constant. The total potential ci> including convection field and corotation 
field is deduced from equations (24) and (1) as; 



(38) 


where k is from equation (28). The potential due to the corotation field is deducec 
from d’j, = /^rcoB(r)d^, because a dipole magnetic field strength at the equator 
is expressed by a radial distance only. In a general distorted magnetic field, the 
corotation field cannot be represented by the gradient of a scalar field. 


The families of equation (37) with various initial conditions c give the particle 
flow patterns, i, e. , particles move along constant total energy lines. Stern 
(1975) studied the characteristic behavior of change in particle energy along the 
dawn-dusk line (0 = ir/2 and 37 t/ 2) using energy conservation law with 7 = 2, and 
gave schematically two kinds of particle flow patterns. His definition of a 
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stagnation point is A (total energy) = 0 which is equivalent to both dx/dt = 0 and 
d(f)/dt - 0. Chen’s calculation (1970) corresponds to the case t = 1 and he also 
used the energy conservation law. 

It is evident from equation (34) that dx/dt = 0 is only satisfied along the dawn- 
dusk line. But, in general the solution for d0/dt = 0 exists at an arbitrary local 
time <p, and henceforth we call this solution a ’’stagnation line", since through 
this line a particle flow direction changes from westward to eastward or vise 
versa, and particle velocities show their local minima along this line. The dawn- 
dusk line is also a stagnation line with respect to x, i. e. , a particle flow direc- 
tion is earthward in the nightside region and is outward in the dayside region. 
Figure 7 depicts one example of flow patterns (solid lines) of particles with M' = 

0. 083 in the orthogonal coordinate system, the ordinate being the normalized 
radial distance x and the abscissa the local time. Thi dashed lines represent 
the stagnation lines and the open arrows indicate the flow directions. The dashed 
areas are trapped regions or forbidden regions inside which any particles with 
the magnetic moment ]U’ never penetrate from the tail regions. As mentioned 
previously, this pattern is independent of electric field strength, only the scale 
factor Rq being a function of electric field strength. In the figure the radius of 
the earth is 0. 2 which corresponds to the stagnation distance Rq = 5a. 

Stagnation Lines and Bast Closed Equlpotentlal Line 


Here we will discuss the stagnation lines for d^/dt = 0 in detail. Equation (31) 
is rewritten (p = sin</>) as; 


f(x) = ~ ^ = px'>' + ^ + 1 - ^ =0. (39) 

CO dt 

And, 

f(x) = (7 + l)p"^ ^ - T (40) 

dt L (t+1)pJ. 

1. For PM’ > 0; 

Since f'(x) has a constant sign, f(x) is a monotonicaEy increasing (p >0 
and m' > 0) or decreasing (p <0 and p' < 0) function, with one solution 
for f(x) = 0. 
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2. For pm’ < 0; 


f • (x) = 0 at X equal to 


and 


Xq = 


^ -6p’ 
d7+ l)Pj 


iJ + 3 


(41) 


f(Xo) 


1p-I v V ^ 6 y 


|p2^4T+ 1)|7+ 3 + 1 


(42) 


2-1. For p >0 and p* <0, f(x) is always positive and there is no solu- 
tion for f(x) = 0. 

2-2. For p <0 and p’ >0, f(x) and f(xQ) are rewritten as 

f(x) = -|p|x'>'+ ^ + 1 -3p7x2, (43) 


and 


where 



7+3 2 


Mo 


7 + 3\ ^ 

7 + iy 



(44) 


(45) 


When 0 < p' < p'q , there always exist two roots (in the case p’ = p’q and 
Ip I =1, two roots become one double root, since f(Xo) = 0). When 
P’o <M’, there is no solution for lp| > Pq and two roots (one double root 
in the case 1 P I = Pg ) for I P 1 ^ Pg, 
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where 


7 + 1 



3. For p = 0; 

one root of X =■>^3^7 independent of 7 . (47) 

4. For p’ = 0; 


1 

one root of X = |p| 7 + 1 when p< 0 , (48) 

and no solution when p > 0 . 

Though the above classifications are complete set of stagnation lines for d0/dt = 
0 , there exists another characteristic magnetic moment P'oo» that is, one equi- 
potential line passes through two stagnation points; one at dusk and another at 
dawn where d^/dt = dx/dt = 0. This magnetic moment m'qq is obtained by a nu- 
merical method from equations (37) and (39). 

Examples of 90° pitch angle particle flow lines and their stagnation lines are cal- 
culated in the case 7 = 2 . Firstly, the last closed equipotential line for p = 0 is 
shown in Figure 8 . This line plays an important role in the formation of the 
plasmapause if we assume the convection electric field to be time-independent, 
because this line is also a boundary of the zero energy particle forbidden region 
inside which particles are trapped and outside which particles take open trajec- 
tories. In general, the equation of this equipotential line is deduced from equa- 
tion (37), c being determined at x = 1 and P = -1. 

px7 + 1 + (7 + l)x - 7 = 0. (49) 

In the case 7 = 1, this is a simple quadratic equation and studied by Chen and 
Wolf (1972). In the case 7 = 2, x is represented in terms of p (= sin<^) as 

(i) For p = 0(<^ = 0 or 7 t), 

X = 2/3. (50) 
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(ii) In the morning side where 0 < p, 


X = 


|(yrT5 + 1)1/3 -(^rTs - 6i 


/3 


In the afternoon (dusk) side where p < 0, 

f- f 6 + 

X = 2 • cos ( 


(51) 


(52) 


and the other solution (x ~>\) is 


/T 


x = 2 • Y 5 cos 


0 + Stt' 


(53) 


where 5 = l/|p| and d = cos ^\/\p \ ; 0 < 0 < tt/2. 


In Figure 8 both cases (7=1 and 7 = 2) are shown, together with Mcllwaln’s in- 
jection boundary by the dashed line (see prev ious sec tion). The dash-dotted lines 
are sta gnation lines discussed above; x = 1/^ \ sin 0 | for 7 = 2 and similarly 

Isin0 I for 7 = 1 when p <0. Note that the shape of Mcflwain’s injection 
boundary normalized by the stagnation point at dusk coincides well with the last 
closed equipotential line with 7 = 2, while the imnormalized distance is far beyond 
that deduced from equations (26) and (32). (See Fig. 6). The ratio of crossing 
points at dusk, midnight and dawn meridians is 

1:2/3: (^2 + 1)’ /^ - (y^- l)i /^ ^ 1 ; 0.67 : 0.60 (54) 

Particle Flow Patterns With 7 = 2 


In the case 7=2, p'oo and are calculated as 

Poo ^ 0-064 

p' = 2“/3 • 5-5/3 _ 0.109. 


(55) 


In order to identify what flow pattern an observed particle takes in the magneto- 
sphere, it is convienient to utilize Figure 9, which illustrates the relation of the 
normalized magnetic moment p' and p with various stagnation distance Lg given 
by equation (36). The upper part of the figure shows a particle energy corres- 
ponding to the magnetic moment indicated in the abscissa. For example in the 
case Lq = 5, p[,q and p[, correspond to magnetic moment p of about 0. 0047 and 
0. 008 [keV/gamma] and their particle energies are about 2. 3 and 3.9 keV at 
L = 4, respectively. The other example is as follows. If the plasmapause 
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position was observed at R = Rp and ^ = 0p , then equation (51) or (52) gives the 
stagnation point Lq at dusk. Using a line labeled , the observed particle 
energy and position give directly the normalized magnetic moment which indi- 
cates which flow pattern that observed particle took. 

A complete set of particle flow patterns classified on the basis on the character- 
istic magnetic moments, i. e. ju' =0, and , is shown in Figures lO(a-g). 
Solid lines represent particle trajectories, flow directions being indicated by 
arrows. Dashed lines are stagnation lines, and A, B, and C are stagnation 
points. The hatched area is a forbidden region and, roughlj?^ speaking, positive 
particles go aroimd eastward where the corotation drift exceeds the gradient-B 
drift and vice versa. The double-hatched area is also a forbidden or trapped 
region, but particles inside this region never make a complete circle aroimd the 
Barth. The size of the earth in the figures is x = 0. 2, i. e. one fifth of the zero 
Bnergy particle stagnation distance Rq at dusk. This can be changed according 
:o the scale factor Rq or Lq as a/Rg = 1 /Lq. 

1. Electrons, M'<0, Figure 10a: 

The flow pattern of electrons is illustrated in Figure 10a for the case 
P* = -0.2. Since directions of gradient-B drift and corotation drift are 
the same sense, i. e. eastward, the radial distance of the stagnation 
point A at dusk is always greater than 1. This indicates that any elec- 
trons injected from the tall or plasma sheet region never penetrate 
inside the plasinapause under the steady state convection electric field. 
The radial distance of the last closed equipotential line that is a boimdarj^ 
of a forbidden region increases as ju’ decreases. 

2. Zero energy particle, M’ = 0, Figure 10b: 

Figure 10b shows the particle flow pattern of zero energy particles. 

The last closed equipotential line (or flow line) has been discussed pre- 
viously in detail (see Fig. 8). Previously it was a misleading thought 
that the boundary inside which the corotation field is dominant and out- 
side which the convection electric field becomes more important is alast 
closed equipotential line or plasinapause location [e. g. Chen, 1970] . 
Strictly the boundary of electric field domin ance is e xpressed in the 
simple form as x = 1 for 7 = 1, and x = {2/\/l + for 7 = 2; the 

latter is represented by die chain circle in the figure, whereas the last 
closed equipotential line is a boundary of the forbidden region (liatched 
area in the figure) given by equations (51) and (52). With respect to the 
p otenti al field the boundary of pot ential field dominance becomes x = 
y/l/\p\ for 7 = 1, and x = ^ 2/1 p I for 7 = 2; the former is identical to 
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the stagnation line in the afternoon to midnight region. The st£^ation 
line indicated by a dashed line is a boundary, in the earth-side region 
of which the eastward drift due to corotation field exceeds the westward 
drift due to convection field for tt < < 2ir, i. e. in the afternoon to mid- 

night. In the region, O<0 <7t, both drifts are eastward. Suffice it to 
\my that this stagnation line is very important for knowing the particle 
flow directions. Since = 0, corotation and convection electric fields 
are simply balanced at the stagnation point A at dusk. Thus tins dis- 
tance Rq represents the strength of the convection electric field. When 
the convection electric field is not steady, that is, time-dependent, a 
last equipotential line is no longer a plasmapause position. 

3. Particles with /x' > 0. 

The flow pattern of positive particle with non-zero magnetic moment 
becomes more complex due to the gradient-B drift, the flow direction 
of wliich is opposite to that of corotation drift. 

(i) 0 < n' < jU^g, Figure 10c; 

The flow pattern of particles with ju' = 0. 029 is illustrated in 
Figure 10c. The size of the forbidden region becomes slightly 
smaller than that for zero energy particles. The stagnation line 
BC appears circled around the earth inside the forbiddesn region 
and its area increases as /x' increases, until the point C coincides 
with the point D. Inside the circle the westward drift mainly due 
to the gradient-B drift dominates the corotation drift, and the par- 
ticle near that stagnation lines does not make a complete circle 
aroimd the earth, just drifting around the stagnation point B. 

(ii) m» = /x'qq. Figure. lOd: 

When jLt' is equal to the characteristic magnetic moment M'qo ~ 

0. 064, the stagnation line inside the forbidden region comes in 
contact with the last closed equipotential line at the point C = D, 
as shown in Figure lOd. The flow pattern is similar to the pre- 
vious one (Fig. 10c), but there is no particle that makes a com- 
plete circle around the earth eastward. The inner equipotential 
line through the point C is an innermost open trajectory of the 
particle which is accessible from the tall region, because the 
forbidden region is separated into two parts by this line when M' 
exceeds M'oo • Note that particles on the eqxiipotential lines 
through A and C take an infinite time to move around the earth. 
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(iii) M{)o < M' <M[) » Figure lOe: 

The flow pattern is shown in Figure lOe for the case jit' = 0. 083. 

The same particle flow pattern is illustrated in Figure 7 in a rec- 
tangular coordinate system to explain the relation of particle flow 
directions with stagnation lines. The separated trapped region 
appears at the dusk side, in which all particles move around the 
inner stagnation point B, In the forbidden region around the earth 
all particles drift westward due to gradient-B drift. Between the 
two separated forbidden regions there is a channel through which 
particles can penetrate closer to the earth along the open trajec- 
tory from the tail. The inmost distance of this trajectory is in 
the dusk meridian and increases as the magnetic moment increases. 

(iv) fJ.' = M'o j Figure lOf: 

When is equal to the characteristic magnetic moment //'q — 

0. 109, the stagnation point A coincides with the point B, i. e. two 
stagnation lines seen in the case 0 <M' < M'q join with each other 
at this point, and the isolated trapped region vanishes away. Fig- 
ure lOf delineates the particle flow pattern of this case. All par- 
ticles inside the forbidden region drift around the earth westward 
because of the dominance of the gradient-B drift over the other 
two drifts. The point A (= B) is an inmost stagnation point at 
dusk; the distance is about x ^^^0.74. 

(V) M'o < M', Figure lOg: 

As the magnetic moment M' increases and exceeds the character- 
istic magnetic moment M'o > there exists only one stagnation point 
at dawn, no stagnation point at dusk, and only one stagnation line, 
as Illustrated in Figure lOg for the case m' = 0.2. The particle 
flow pattern for n' > (jl'q is apparently similar to those of a zero 
energy particle or electrons, if the directions of dawn and dusk 
are exchanged around the noon-midnight meridian. The stagna- 
tion line, however, is distinctly different for each case; that is, 
there always exists particles which change their flow directions 
four times in the case ju’ > M'o » whereas there are no such elec- 
trons or zero energy particles. In other words, the stagnation 
line for M' ^ in the afternoon to midnight region never vanishes, 
while there is ho stagnation line for M’ ^ 0 in the midnight to after- 
noon region. 
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VI. TIME-DEPENDENT PARTICLE FLOWS 


The observational data of energetic particles associated with substorms show 
complicated but clear particle enhancements around or inside the plasmapause 
[Smith and Hoffman, 1974] . Trajectories of those particles which are presum- 
ably injected from the plasma sheet or tail region and drift in the magnetosphere 
are strongly dependent on particle energy and pitch angle due to a combination of 
convection, corotation, and gradient-curvature drifts. 

To understand the fundamental particle motion in the magnetosphere, the 90° 
pitch angle particle motion is studied in detail, and a complete set of flow pat- 
terns is presented, taldng into account the stagnation lines along which d<^/dt = 

0; the stagnation line for dx/dt = 0 is a dawn-dusk meridian line. Although those 
patterns explain the energy-dependent and complicated trajectories of particles, 
the drift time along the trajectory can not be deduced from the energy conserva- 
tion law and is not indicated in the figures. To interpret the observed enhance- 
ments of injected particles, more important is a dynamical motion of the particles 
or a particle tracing with respect to the time-phase of a substorm and the obser- 
vational time. 

Using the three simultaneous differential equations in (C-1) to (C-4), particle 
trajectories in the magnetosphere are simulated. In Figures ll(a-h) are illus- 
trated some examples of results for 90° pitch angle particles, which display the 
energy dependence on the particle trajectories. The initial particle source lo- 
cation is assumed to be at L = 10 and 20:00 MLT < 0 < 4:00 MLT at 0. 1 hour 
intervals. The initial particle energies range from -0. 4keV to 51. 2keV for 
Figures ll(a-l). The strength of the convection electric field is represented by 
the stagnation distance Lq = 5, which corresponds to A = 0. 057 [ mV/m* R^] and 
an electric field intensity of 0. 57 [mV/m] at L = 10 in the midnight. From 
equation (26) this is the case for Kp — 4. 1. The particle positions are plotted at 
0. 1 hour steps by dots and at 1 hour steps by asterisks, within 20 hours and 
L < 10. These particle traces can be easily classified into the categories of par- 
ticle flow patterns given in Figure 10, by means of the relation of ju' (normalized 
magnetic moment) and W (particle energy) or p in Figure 9; (a) to Figure 10(a), 

(b) to Figure 10(b), (c) to Figure 10(c), (d) to Figure 10(e), and (e)^/ (1) to Fig- 
ure 10(g). In these particle trajectories, therefore. Figure 11(d) shows the in- 
most accessible trajectory to the earth if the infinite time of particle travelling 
is considered. But, as is evidently understood in these fignres the inmost acces- 
sible distance strongly depends on particle travelling time from the source region 
as well as particle energy. This is the main cause of the formation of a 'nose' 
structure inside the plasmapause observed by -A [Smith and Hoffman, 1974]. 


20 



In the midnight region, particles with energies less than about lOkeV at L = 10 
come from L = 10 to L = 9 in about 6 min, to L = 8 in about 18 min, and L = 7 
in about 36 min, which are determined mainly by the convection drift velocity. 

A rough estimate of the convection speed 0.1 Rj^/min beyond the synchronous 
orbit from ATS-1 magnetic field variations derived by Parks et al. (1971) are 
in the range of this result: 0. 17 R^/min to 0. 06 Rp/min. Also, the electric 
field intensity at synchronous orbit used in this calculation is 0. 376mV/mwMch 
coincides well with the value of 0.36 mV/m by Shelley et al. Note that this con- 
vection electric field is a function of magnetic activity and position given in Equa- 
tions (24) and (25). 

Ambiguity in the source location of the particles does not affect much die particle 
trajectories given in the figures since the front of the particle trajectories, i. e. 
all particle locations at a fixed time beyond the synchronous orbit, is almost ir- 
respective of particle energy, pitch angle and charge. Therefore, only the time 
marks need be shifted so that the time of the particles at a source location is a 
starting time. For example in the case of a particle source location at L = 8 at 
midnight, the first asterisk shows the location of a particle at the time of +0. 7 
hour, etc. 

Although these trajectories are immutable with respect to the convection electric 
field strength as previously mentioned, the actual particle locations are very 
sensitive to the change in distance of the stagnation point, because of changes of 
both scale factor and trajectory due to change of magnetic moment (see equations 
(32) and (36)). For the higher energy particles shown in Figure ll(k-l), the tra- 
jectories look very different. There are no particles at the dawn meridian, and 
futhermore particles are unable to come inside the synchronous orbit, but are 
swept out towards the dayside magnetopause very rapidly. 

Fignres 12(a-d) illustrate the difference of particle trajectory due to pitch angles. 
The particle initial energy at L = 10 is 0. 566 keV and their pitch angles are (a) 

0°, (b) 30°, (c) 60°, and (d) 90°. As is shown in Figure 3 the energies of particles 
with different pitch angles change with L at different energization rates, and par- 
ticles take quite different trajectories especially inside the synchronous orbit. 
This means that particles with the same energy but different pitch angles starting 
from the same source location are dispersed to different locations according to 
their energy and pitch angle. This may be the main cause of pitch angle aniso- 
tropies found in the vicinity of the proton 'nose' structures associated with sub- 
storms [Ejiri et al, , 1975], For most of substorm associated particles with 
energies of about IkeV to 30keV inside the plasmapause in the dusk side, parti- 
cles with large pitch angles penetrate to lower L-values than particles with small 
pitch angles. This is clearly seen in Figure 13 which shows the locations of par- 
ticle fronts for particles with 0° and 90° pitch angles; other pitch angle particle 
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fronts are located between these two lines. From this figure we can speculate 
a quite different pitch angle anisotropy in die morning side, i. e. inverse .aniso- 
tropy to that in the evening side. The distance between different pitch angle par- 
ticle fronts changes with travelling time and depends on the local time meridian 
to be considered as well. 

The convection electric field used in this calculation is time -independent. Re- 
cently Roederer and Hones (1974) developed a time-dependent electric field model 
to interpret die ATS- 5 particle data. Though the electric field as a function of 
Kp is time-dependent, the time resolution of 3 hours is not sufficient to interpret 
the particle enhancements during substorm expansion phase. 

The detailed descriptions of observational data of Explorer 45 (S^ -A) and their 
interpretations using particle tracings developed in this paper will be presented 
in a separate paper. 


VII. ACKNOWLEDGMENTS 

I wish to express my sincere thanks to Drs. R. A. Hoffman and P, H. Smith 
for their encouragement and valuable discussions, especially for providing un- 
published Explorer 45 particle data that motivated much of this work, and to 
Drs. D. P. Stern, A. J. Chen, M. Sugiura, and K. Maeda for their comments 
and suggestions. I am grateful to Christine Gloeckler, Christine J. McQuillan 
and Dr. Narindra K. Bewtra for assistance in the computer programming and 
displays. I am also indebted to NAS/NRC for sponsoring my Post- Doctor 
Research Associateship in the Laboratory for Planetary Atmospheres, NASA/ 
GSFC. 


22 



VII. REFERENCES 


Cauffmati, D. P. and N. C. Majmard, "A Model of the Effect of the Satellite 
Photosheath on a Double Floating Probe System", J. Geophys. Res. , 79 , 
2427-2438, 1974. 

Chen, A. J. , "Penetration of Low-energy Protons Deep into the Magnetosphere", 
J. Geophys. Res., 2458-2467, 1970. 

Chen, A. J. , and R. A. Wolf, "Effects on the Plasmasphere of a Time-varying 
Convection Electric Field", Planet. Space Sci. , 483-509, 1972. 

Chen, A. J. , and J. M. Grebowsky, "Plasma Tail Interpretations of Pronounced 
Detached Plasma Regions Measured by OGO 5", J. Geophys. Res., Td, 
3851-3855, 1974. 

Chen A. J. , and J. M. Grebowsky, and H. A. Taylor, Jr., "Dynamics of Mid- 
latitude Light Ion Through and Plasma Tails", J. Geophys. Res. , 

968-976, 1975, 

Chen, A. J. , and D. P. Stern, "Adiabatic Hamiltonian of Charged Particle Motion 
in a Dipole Field", J. Geophys. Res., 80, 690-693, 1975. 

Cowley, S. W. H. , and M. Ashour-Abdalla, "Adiabatic Plasma Convection in a 
Dipole Field; Variation of Plasma Bulk Parameters with L", Planet. Space 
Sci. , 1527-1549, 1975. 

DeForest, S. E. , and C. E. Mcllwain, "Plasma Clouds in the Magnetosphere", 

J. Geophys. Res. , 3587-3611, 1971. 

Ejiri, M. , R. A. Hoffman, and P. H. Smith, "Model Calculations of Particle 
Traces in the Magnetosphere - the Nose Structure of Ring Current Proton 
Enhancements Observed by Explorer 45 (S^ -A)", Trans. Am. Geophys. 

Union, EOS, No. 12, 1047, 1975. 

Frank, L. A., "Plasma in the Earth’s Polar Magnetosphere", J. Geophys. Res., 
76, 5202-5219, 1971. 

Grebowsky, J. M. , Y. Tvilunay, and A. J. Chen, "Temporal Variations in the 
Dawn and Dusk Midlatitude Through and Plasmapause Positions", Planet. 
Space Sci. , 1089-1099, 1974. 

Gurnett, D. A. , and L. A. Frank, "Observed Relationships Between Electric 
Fields and Auroral Particle Precipitation", J. Geophys. Res., 

145-170, 1973. 


23 



REFERENCES (continued) 

Heikkila, W. J. , "Outline of Magnetospheric Theory", J. Geophys. Res., 7^, 
2496-2500, 1974. 

Hamlin, D. A. , R. Karplus, R. C. Vik, and K. M. Watson, "Mirror and 
Azimuthal Drift Frequencies for Geomagnetic ally Trapped Particles", 

J. Geophys. Res. , 1-4, 1961. 

Heppner, J. P. , "Electric Field Variations During Substorms; OGO 6 Measure- 
ments", Planet. Space Sci. , 1475-1498, 1972. 

Heppner, J. P. , "High-latitude Electric Fields and the Modulations Related to 
Interplanetary Magnetic Field Parameters", Radio Sci. , 8, 933-948, 1973. 

Hilton, H. H. , "L Parameters, a New Approximation", J. Geophys. Res. , 7£, 
6952-6954, 1971. 

Kavanagh, L. D. Jr. , J. W. Freeman, Jr. , and A. J. Chen, "Plasma Flow in 
the Magnetosphere", J. Geophys. Res. , 73, 5511-5519, 1968. 

Kivelson, M. G. , and D. J. Southwood, "Local Time Variations of Particle 

Flux Produced by an Electrostatic Field in the Magnetosphere", J. Geophys. 
Res., 80, 56-65, 1975. 

Konradi, A. , C. L. Semar, and T. A. Fritz, "Substorm-injected Protons and 
Electrons and the Injection Boundary Model", J. Geophys. Res. , M, 

543-552, 1975. 

Lenchek, A. M. , S. F. Singer, and R. C. Wentworth, "Geomagnetic ally Trapped 
Electrons from Cosmic Ray Albedo Neutrons", J. Geophys. Res. , 66 , 
4027-4046, 1961. 

Lew, J. S. , "Drift Rate in a Dipole Field", J. Geophys. Res., 66, 2681-2685, 
1961. 

Lezniak, T. W. , and J. R. Winclder, "Experimental Study of Magnetospheric 
Motions and the Acceleration of Energetic Electrons during Substorms", 

J. Geophys. Res., 75, 7075-7098, 1970. 

Longanecker, G. W. , and R. A. Hoffman, "S^ -A ^acecraft and Experimental 
De:-;cription", J. Geophys. Res., 78, 4711-4717, 1973. 


24 



REFERENCES (continued) 


Mauk, B. H. , and C. E. McRwain, "Correlation of Kp with the Substorm-injected 
Plasma Boimdary", J. Geophys. Res. , 3193-3196, 1974. 

Maynard, N. C. , and D. P, Cauffman, "Double Floating Probe Measurements 
on S^ -A", J. Geophys. Res. , 4745-4750, 1973. 

Maynard, N. C. , and A. J. Chen, "Isolated Cold Plasma Regions: Observations 
and Their Relation to Possible Production Mechanisms", J. Geophys. Res. , 
1009-1013, 1975. 

McRwain, C. E. , "Plasma Convection in the Vicinity of Geosynchronous Orbit", 
Earth's Magnetospheric Processes, ed. B. M. McCormac, 268-279, 

D. Reidel Pub. Co. Dordrecht-HoUand, 1972. 

McRwain, C. E., "Substorm Injection Boundaries", paper presented at the^ 
Chapman Memorial Symposium, NOAA, Boulder, Colorado, June, 1973. 

Nakada, M. P. , J. W. Dungey, and W. N. Hess, "On the Origin of Outer -belt 
Protons, 1", J. Geophys. Res., 7^, 3529-3532, 1965. 

Nishida, A. , "Formation of Plasmapause, or Magnetospheric Plasma Knee, by 
the Combined Action of Magnetospheric Convection and Plasma Escape from 
the Tsiil", J. Geophys. Res. , n, 5669-5679, 1966. 

Northrop, T. G. , and E. Teller, "Stability of Adiabatic Motion of Charged Par- 
ticles in the Earth's Field", Phys. Rev. 117 , 215-225, 1960. 

Northrop, T. G. , The Adiabatic Motion of Charged Particles , New York, Inter- 
science Publishers, 1963. 

Parks, G. K. , E. H. Gardner, C. S. Lin, R. E. NeweR, and P. J. Wehrenberg, 
"Interplanetary Control of Magnetosphere Dynamics", Correlated Interplane- 
tary and Magnetospheric Observations, ed. D. E. Page, 49-60, D. Reidel 
Pub. Co. Dordrecht-Holland, 1974. 

Roederer, J. G. , Dynamics of GeomagneRcaRy Trapped Radiation , Physics and 
Chemistry in ^ace 2, New York, ^ringer- Verlag, 1970. 

Roederer, J. G. , and E. W. Hones, Jr. , "Motion of Magnetospheric Particle 
Clouds in a Time-dependent Electric Field Model", J. Geophys. Res. , 
1432-1438, 1974. 


25 




REFEKENCES (continued) 

Schield, M. A. , and L. A. Frank, "Electron Observations Between the Inner 
Edge of the Plasnaa Sheet and the Plasrnasphere", J. Geophys. Res. 7 ^, 
5401-5414, 1970. 

Schulz, M. , "Approximate Second Invariant for a Dipole Field", J. Geophys. 
Res., 3144-3148, 1971. 

Shelley, E. G. , R. G. Johnson, and R. D. Sharp, "Plasma Sheet Convection 
Velocities Inferred from Electron Flux Measurements at Synchronous 
Altitude", Radio Sci. , No. 2, 305-313, 1971. 

Smith, P. H. , and R. A. Hoffman, "Ring Current Particle Distributions during 
the Magnetic Storms of December 16-18, 1971", J. Geophys. Res. , 7^, 
4731-4737, 1973. 

4» 

Smith, P. H. , and R. A. Hoffman, "Direct Observations in the Dusk Hours of 
the Characteristics of the Storm Time Ring Current Particles During the 
Beginning of Magnetic Storms", J. Geophys. Res. , 966-971, 1974. 

Southwood, D. J. , and M. G. Kivelson, "An Approximate Analytic Description 
of Plasma Bulk Parameters and Pitch Angie Anisotropy Under Adiabatic 
Flow in a Dipolar Magneto spheric Field", J. Geophys. Res. , 2069- 

2073, 1975. 

Stern, D. P. , "Models of the Earth’s Electric Field", NASA/GSFC, preprint, 
X-602-74-159, May, 1974. 

Stern, D. P. , "The Motion of a Proton in the Equatorial Magnetosphere", 

J, Geophys. Res., 89, 595-599, 1975. 

Sugiura, M. , and S. Chapman, "Tables and Figures Related to the Geomagnetic 
Dipole Field", NASA/GSFC, document X-612-66-546, December, 1966, 

Taylor, H. E. , and E. W. Hones, Jr. , "Adiabatic Motion of Auroral Particles 
in a Model of the Electric and Magnetic Fields Surrounding the Earth", 

J. Geophys. Res., 70, 3605-3628, 1965. 

Vasyliunas, V. M. , "A Survey of Low-energy Electrons in the Evening Sectors 
of the Magnetosphere With OGO 1 and OGO 3", J. Geophys. Res. , 
2839-2884, 1968. 


26 



REFERENCES (continued) 


Volland, H. , "A semiempiric al Model of Large-Scale Magneto spheric Electric 
Fields”, J. Geophys. Res., 171-180, 1973. 

Williams, D. J. , R. A. Hoffman, and G. W. Longanecker, "The Small Scientific 
Satellite (S^) Program audits First Payload”, NASA/GSFC, Preprint 
X-612-68-418, 1968. 


27 



APPENDIX A 


The following symbols are used in the text: 

a earth’s radius (6.371 x 10^ m) 

equatorial pitch angle 

s' dipole magnetic field; B = ky/L^ at equator 

Bm B at X = 

E electric field; E =-V^^[V/m] 

J second (longitudinal) adiabatic invariant 

Icq magnetic field strength atL = l (0.31x Wb/m^) 

L Mcllwain’s L- value; E = aL 

X latitude 

Xj^^ mirror latitude 

M first adiabatic invariant (magnetic moment); fx = p^/2m 

p particle momentum; p = mv 

tl> electric potential 

(P local time in radian 

q, m particle charge (negative for electrons) and mass 

E particle guiding center position at equator 

E(, stagnation distance at dusk for zero energy particle; Eq = uLq 

s length along the magnetic field line 

Sj^ , s^ mirror points 

t time 

V particle velocity 

W particle lane tic energy 

u) angular velocity of the earth 

(i and II denote the components perpendicular and parallel to the magnetic 
field line , respectively) 

NOffi.wa 



APPENDIX B 


To obtain an accurate value of integration of equation (6), the integrand of which 
has a singularity, the integral f(aQ) is divided into two parts, i.e. 

- H 

F(z)dz + 

0 


/' 

•' 1 - H 


F(z)dz 


(B-1) 


where 


z = sin X/sin . 


The integrand of the second term is expanded around tlie singular point as 


F(z) = g,(a) 



+ §2 


(a)h''^^ + 0(h) 



(B-2) 


where 


z = 1 - h;he[l -H, 1], 
a = sinXj^(ajfc 1), 

g,(a) = (1 +3a-)yi - a V\/3(3 - 5a2) , 

and 

g 2 (a) = (60a^ + 19a^ - 12a^ + 3)/4(l +3a^)(l -a2)(3 + 5a“). 

Then, the integration of F(z) in a region [1-H,1] becomes 

^i(%) = f F(z)dz = gj(a)|2H’/^ + 2/3g2(a)H3/2 + (B-3) 

The integral f^ (ce^) can be computed with a specified accuracy by usual numerical 
integration methods, since the functional values of F(z) and its n~th derivative 
have finite values in the region [0, 1-Hl. From equation (B-3), we have the 
value of f(Q!o) with an accuracy of ; for instance if H is 10"^ , then ~ 

10"'^ It is further recommended to divide the region [0, 1-H] into several 
regions and to use a sufficiently small step corresponding to each integral region. 


Preceding 


Rage 
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—8 

keeping the same accuracy. In order to get the accuracy of 10 by the Hermitian 
fourth order integration method, the region is divided as [1 -10"'^ , 1 -10“^'^ ^ ^ ; 

n = 0, 1, . . . 5, and a numerical integral step of Ah = 0. O/IO’"* ^ is used. This 

method reduces the nurhber of calculating steps to the amount 10"^, compared 
with the case of a constant integral step used throughout the region. 

The integral I(<^o) singularity in the integrand but its derivative has the 

same kind of singularity as F(z). The above method therefore, should be used. 

Tins method can be applicable to any integral in the case of its integrand or its 
n-th derivative having any singular points. Sigiura and Chapman (1966) developed 
another method of integration and gave the functional values of and 1 ( 0 ;^) with 
an accuracy of 10~^ at 1° intervals of mirro r latitude. 
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APPENDIX C 


For practical use, equations (30), (31), and (21) are transformed, with change 
of the units (y?^0) into 


dt 



xT + 2 

COS0', 

7 


dd> . U' 

= 1 ^ sin0' + C2G(y) 


dt 




and 


_ _L . Z . ^ 

dt “ T ’ 7 ’ ’ dF’ 


(C-1) 


(C-2) 


(C-3) 


where dx/dtin l/[hour:UT], d0/dt in [hour:LT] / [hour:UT], cp' = Cj0, U’ = 
M/q*y^ in [kV/gamma], and C, = -tt/12 and C 2 = -3 x 10^ ^ /a.^ oj — -lOlQ. 


In the case y = 0, from equation (16) the third term of equation (C-2) becomes 


C W. 




(C-4) 


where Wj is an initial energy in keV, Xj an initial position at t = 0, and 
C3 = -2 X loVa^coko ci -0. 022. 


The above simultaneous differential equations are solved numerically with 
Hamming’s modified predictor-corrector method, an initial condition being (x, 
y> W; t= 0) = (X|, 0j, yj, Wj). Using approximate formulas for f(y), I(y), 
and G(y) in equations (17), (20), and (22), the solution (x, 0, y; t) is obtained with 
sufficient accuracy of 5 x 10"^ , the time step being automatically changed. TMs 
accuracy is determined from the observational data by Explorer 45 (S ^ -A) which 
show that the difference in position (x, 0) of the particle enliancement between dif- 
ferent pitch angle particles is about 10“^ (see Fig. 1). The resolution of calcu- 
lated particle position, therefore, is required to be about 10"^. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2a. 


Figure 2b. 


Figure 3. 


Figure 4a - c. 


Figure 5a. 


Figure 5b. 


Figure 6. 




Proton spectrograms of orbit 281 observed by Explorer 45 
(S' -A) on Feb. 13, 1972. The plasmapause positions during 
the outbound and inbound portions of the satellite orbit were iden- 
tified by the DC electric field experiment [Maynard and Cauffman, 
1973; Cauffman and Maynard, 1974] at L = 5. 1 and L = 4. 8, re- 
spectively. The satellite position was about L = 5 at 19:20 UT out- 
bound and at 21:35 UT inbound. The magnetic local time at apo- 
gee was 19:10 MLT. The typical nose structure in each pitch 
angle spectrogram is seen inside the plasmapause. 

Change in equatorial pitch angle with L-value. The equatorial 
pitch angles indicated on the right side are values at L = 10. 

Variations of equatorial pitch angles at various distances indi- 
cated by L-values as parameters. The particle which changes 
its pitch angle beyond the loss cone limit falls into the loss cone. 

Changes in particle energy due to radial drift motion. Indicated 
parameters are initial pitch angles at L = 10. A pitch angle also 
changes with L as seen in Figure 2, except at 90° and 0°. 

The functional values of f(«o)> I(^q) ^vid GCceg) as a function of 
equatorial pitch angle or colatitude of mirror point: f(0°) = 
0.1380172998, 1(0°) --- 0.2760345996, C(0°) = 2/3, f(90°) = 
0.7404804896, 1(90°) = 0, and G(90°) = 1. 

Relative accuracies of the functions f(y) (dashed line) and I(y) 
(solid line), with mirror colatitude computed from equa- 
tions (17) and (20). The results from Schulz's approximations 
are also shown. 

Relative accuracy of the fimction G(y), with mirror colatitude 
computed from equation (22). The result from SclixiLz's 
approximation is also shown. 

Kp-dependence of plasmapause position in the case of time- 
independent electric field model with 7 = 2. is the dis- 

tance of the stagnation point at dusk in equation (33), and 
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FIGUKE CAPTIONS (continued) 


Figure 6. 
(continued) 


Figure 7. 


Fig-ure 8. 


Figure 9. 


Figure 10a - 


Figure 11a - 


Lj^,j [ j) i(; j., j is the distance of crossing point of the last closed 
equipotential line at the midnight meridian from equation (54). 

represents the radial distance of Mcllwain's injection boun- 
dary at each local time; = (122 -10 Kp) / (0 -7. 3) where Rt, 
is L -value and 0 the local time in hours. 

Example of particle trajeetories (solid line) in the case M' = 

0. 083. Dashed lines represent the stagnation lines through 
which a particle changes its drift direction. The horizontal 
open arrows indicate directions of azimuth drift and the vertical 
open arrows indicate directions of radial drift. 

Last closed equipotential lines of zero energy partiele with 
7 = 1 and 2. The dashed line from dusk to midnight shows 
Mcllwain's injection boundary normalized by the stagnation 
distance at dusk. The dash-dotted lines are stagnation lines. 

The relation of normalized magnetic momeiat/x' and |u[ eV/gammal 
with various stagnation distance Lq as parameters. Two charac- 
teristic normalized magnetic moment ^ 0. 109 and 
0. 064 are indicated by dashed lines. The particle having a mag- 
netic moment between these lines (indicated by the dark solid 
line) takes a flow pattern depicted in Figure 7 or Figure lOe. 

The upper part shows the particle energy as a function of radial 
distance L, corresponding to the magnetic moment indicated in 
the abscissa. 

Particle flow patterns (a) electrons, g' = -0.2 <0, (b) zero 
energy particle, ja' = 0, (c) 0 <M' = 0. 029<M'oo » (d) 

0.064, (e) <m’ = 0. 083 <m'q, (f) M' = M'o —0.109, aiid 

(g) M’o <M' = 0.2. 

The 90° pitch angle particle tracing computed from equations 
(C-1) to (C-4). Initial particle positions are at L = 10 and 
20:00 MLT < 0 <4:00 MLT at 0. 1 hour interval. Particle posi- 
tions are represented by dots at 0. 1 hour (= 6 min) steps and by 
asterisks at 1 hour steps witMn 20 hours and L <10. Stagnation 
distance is Lq = 5, constant. The initial particle energies Wj 
at L = 10 are; (a) -0.4keV (electron), (b) 0, (c) O.lkeV, 

(d) 0. 2keV, (e) 0.4keV, (f) 0. 8keV, (g) 1.6keV, (h) 3.2keV, 

(i) 6. 4keV, (j) 12. 8keV, (lc)25.6keV, and (L) 51. 2keV. 
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FIGURE CAPTIONS (continued) 


Figure 12a- 


Figure 13. 


. The trajectories of particles with the energy of 0. 566 keV at 
L = 10, and pitch angles of (a) 0°, (b) 30°, (c) 60°, and (d) 90°. 
See the caption of Figure 11. 

The difference of particle flow fronts due to pitch tmgies, 0° 
and 90°. Initial particle positions and electric field intensity 
are the same as Figure 12. Particle energies at L = 10 are 
1. 6 keV for 0° particles and 0. 8keV for 90° particles, and tliese 
particle energies are approximately equal at L = 3. 5. 
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Figure 1 (Continued). 
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Figure 6. 
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(b) W, = OkeV (M*0) 
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•Figure 11. 
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(g) W,=16 keV (^ = 52 eV 7) 
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Figure 11 (Continued). 
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(a) Oo* 0“ ( W, * 0 566 KeV./i = 0) 



12 

(b)ao*30°(W,=0566 KeV>i = 46eV/y) 



(c)ao*60“(W| = 0666KeV/i. = I37eV /) 



12 

(d)a„=90 (W, = 0 566 KeV /i= 182 eV Y) 


Figure 12 



1.6 keV 16.9 keV 0° 
0.8 keV 18.6 keV 90° 
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